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Abstract 

We use the Fuchsian algorithm to study the behavior near the sin- 
gularity of certain families of C/(l) Symmetric solutions of the vacuum 
Einstein equations (with the U{1) isometry group acting spatially). 
We consider an analytic family of polarized solutions with the max- 
imum number of arbitrary functions consistent with the polarization 
condition (one of the "gravitational degrees of freedom" is turned off) 
and show that all members of this family are asymptotically velocity 
term dominated (AVTD) as one approaches the singularity. We show 
that the same AVTD behavior holds for a family of "half polarized" 
solutions, which is defined by adding one extra arbitrary function to 
those characterizing the polarized solutions. (The full set of non- 
polarized solutions involves two extra arbitrary functions). We begin 
to address the issue of whether AVTD behavior is independent of the 
choice of time foliation by showing that indeed AVTD behavior is 
seen for a wide class of choices of harmonic time in the polarized and 
half-polarized {U (f ) Symmetric vacuum) solutions discussed here. 



1 



1 Introduction 



During the last few years, our understanding of the behavior of cosmological 
solutions^] near their big bang singularities has increased significantly. There 
is more and more evidence, both numerical [0] @] 0] [0 and analytical 
HI m [0 [Q ' ^h^t wider and wider families of such spacetimes exhibit either 
asymptotically velocity term dominated ("AVTD") behavior, or Mixmaster 
behavior, in a neighborhood of their singularities. 

A solution has AVTD behavior if the metric (when expressed in suitable 
gauge) asymptotically behaves at each spatial point like a Kasner spacetime 
metric, with the Kasner parameters generally varying from point to pointQ. 
A solution exhibits Mixmaster behavior if, instead of asymptotic Kasner 
evolution at each point, one asymptotically has at each point the oscillatory 
behavior of bouncing from one Kasner epoch to another, as seen in Bianchi 
IX cosmologies. The idea that a generic solution should have one or the 
other of these behaviors was first suggested by Belinskii, Khalatnikov, and 



Lifschitz [jT2|, |T3[, and is therefore known as the "BKL conjecture." 

For spacetimes with a U{1) isometry group and spatial topology, there 
is strong numerical support for the BKL conjecture. Extensive numerical 
studies of these spacetimes indicate that if the Killing field generating the 
isometry is hypersurface orthogonal (these are the so-called "polarized so- 
lutions," since one of the two gravitational degrees of freedom is essentially 
turned off), then AVTD behavior is seen 0; while for generic U{1) Symmet- 
ric spacetimes with spatial topology, one finds Mixmaster behavior 0. 
These results agree with the conclusions of earher studies of U{1) Symmetric 



spacetimes [|T4| , which used formal expansions and heuristic analysis to argue 
that the BKL conjecture likely holds for them. 

In this work, we rigorously prove some of these results, using Fuchsian 
methods. More specifically, we show that for rather general sets of ana- 
lytic data for polarized U{1) Symmetric solutions, the solutions are AVTD. 
Furthermore, we identify a class of U{1) Symmetric solutions, intermedi- 
ate between the polarized and the generic classes, which is AVTD as well. 
The solutions in this class, which we call "half-polarized," are intermediate 
in the sense that, modulo gauge freedom and the constraint equations, the 



^Here we follow Bartnik's convention Q of calling a spacetime a "cosmological solution" 
if, in addition to being a solution of Einstein's equations, it is globally hyperbolic, with 
compact Cauchy surfaces. 

^See Section 3 for a more precise definition of AVTD behavior. 
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polarized solutions are characterized by two arbitrary functions on T^, the 
half-polarized by three such functions, and the nonpolarized by four. So in 
a rough sense, what we find is a distinct way to turn on half of the second 
degree of freedom at each point, without removing the AVTD behavior. We 
further show that the class of U{1) Symmetric solutions with AVTD behavior 
includes spacetimes which are neither polarized nor half-polarized, since we 
can find such solutions by applying SL{2, R) parametrized Geroch transfor- 
mations to the half-polarized solutions and still retain AVTD behavior. 

Fuchsian methods have become an important tool in recent years for 
proving the existence of AVTD behavior in cosmological solutions of Ein- 
stein's equations. Kichenassamy and Rendall first used these techniques, in 
1998, to establish AVTD behavior in (unpolarized) Gowdy spacetimes [^. 
Subsequent work |^ showed that in polarized Symmetric solutions with 
nonvanishing twist (the vanishing of the twist characterizes the Gowdy solu- 
tions as a subfamily of the Symmetric cosmological spacetimes), AVTD 
behavior is found. Especially interesting is the very recent result JlTI which 
uses Fuchsian methods to show that AVTD behavior occurs in cosmological 
solutions (spatially T^) of the Einstein-scalar field equations with no assumed 
symmetries. 

In none of these works, including the work discussed here, do the Fuchsian 
methods show that all solutions in the family under consideration (general 
Gowdy, polarized Symmetric, polarized and half-polarized U{1) Symmet- 
ric, etc) necessarily have AVTD behavior. Rather, the Fuchsian methods 
show that in each family of solutions, there is an analytic subfamily which is 
defined by the same number of free functions as the full family, and whose 
members all exhibit AVTD behavior. Thus, in a rough sense, Fuchsian meth- 
ods show that AVTD behavior is stable, occurring in "open subsets" of the 
given full family. It is expected, based on numerical evidence, that AVTD 
behavior is generic in these families; however, this has not been proven yet 
in any family except for the polarized Gowdy solutions. 

The polarized and half-polarized U{1) Symmetric spacetimes with which 
we work here comprise the least restricted family of vacuum solutions in 
which AVTD behavior has been proven to exist. It is useful to note that the 
recent work showing that AVTD behavior occurs in Einstein-scalar field so- 
lutions with no symmetry does not imply our U{1) results. Indeed heuristic 
arguments first presented by Belinskii, Khalatnikov, and Lifschitz [|13| indi- 
cate that the presence of a nonvanishing scalar field is likely to remove the 
inevitability of the "potential bounces" which lead to Mixmaster behavior. 
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with AVTD behavior resulting. It is expected that vacuum solutions with 
no isometrics are generically not AVTD. 

It would, of course, be very nice if Fuchsian type methods could be used to 
rigorously prove that Mixmaster behavior occurs in those families of solutions 
(like the magnetic Gowdy spacetimes [Q, the nonpolarized Symmetric 
spacetimes 0, and nonpolarized U{1) Symmetric spacetimes 0) in which 
Mixmaster behavior has been observed numerically. However, it is not known 
how to do this. Indeed, only very recently has Mixmaster behavior been 
rigorously verified in spatially homogeneous spacetimes [|l| 



In all Fuchsian studies of cosmological solutions prior to ours, the choice 
of spacetime foliation has been more or less rigidly determined by the family 
of spacetimes and the analysis. For example, for the Gowdy and for the 
Symmetric spacetimes, the analysis of solutions is simplest if one uses the 
"areal" (or "Gowdy") foliation ||T^ |]TB[, in which the t = constant hypersur- 
faces each consist of orbits of constant area; this foliation is (up to scaling) 
unique. In our present study of polarized U{1) Symmetric spacetimes, how- 
ever, the analysis is carried out using a harmonic time foliation, as detailed 
below. Harmonic time is not unique, and so we have the opportunity here to 
consider the issue of whether, and in what sense, observed AVTD behavior 
depends on the choice of foliation. We show that if, in a fixed U{1) Symmet- 
ric solution (T^ X R,g), AVTD behavior is seen using one choice of harmonic 
time, then there is a full (two free functions on T^) family of other choices of 
harmonic time for (T^ x R,g) such that AVTD behavior is seen using each of 
them as well. We also discuss the extent to which we expect AVTD behavior 
to be seen using other foliations. 

The bulk of this paper is devoted to the verification, using Fuchsian meth- 
ods, that AVTD behavior occurs in polarized and half-polarized U{1) Sym- 
metric vacuum solutions on x R. We begin by reviewing (in Section 2), 
the form of the metric and the form of the field equations we shall work with 
here. We obtain these forms by starting with the generic metric for U{1) 
Symmetric spacetimes, and then imposing certain gauge conditions, and cer- 
tain restrictions on the fields which are preserved by the field equations. We 
write the field equations in canonical, Hamiltonian form. We note in Section 
2 that in setting one of the fields together with its conjugate momentum 
to zero, we define the polarized solutions. We further note that the half- 
polarized solutions cannot be characterized by restrictions on the fields at 
an arbitrary finite time; rather, they are characterized by certain asymptotic 
conditions near the singularity. 
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In Section 3, we write the asymptotic (near the singularity) ansatz for the 
polarized solutions, and then for the half-polarized solutions. The asymp- 
totic ansatz for the polarized fields (forgetting the constraints) involves eight 
free functions (the "asymptotic data") on together with eight remainder 
functions depending on time as well as on T^. (There are eight rather than 
two because the constraint equations have not been imposed, and some gauge 
freedom remains.) Substituting this ansatz into the Einstein evolution equa- 
tions, one obtains evolution equations for these remainder functions (with 
the free functions on appearing as parameters). We show in Section 3 
that for all choices of the asymptotic data which satisfy certain linear in- 
equalities, these evolution equations take the special Fuchsian form. It then 
follows from a result of Kichenassamy and Rendall that for any such 
choice of the asymptotic data, there is a unique solution for the remainder 
functions, decaying to zero near the singularity. This behavior guarantees 
that solutions in the (polarized) ansatz form are AVTD. The analysis for so- 
lutions in the half-polarized asymptotic ansatz form is similar, with similar 
conclusions. 

The discussion in Section 3 ignores the constraint equations. In Section 
4, we show that if the asymptotic data in the asymptotic ansatz expressions 
satisfy certain asymptotic constraints (|67D — (|^), then the corresponding 
solution satisfies the Einstein constraints. The converse is true as well; so 
the Einstein constraints are (at least among solutions of the evolution equa- 
tions) equivalent to these asymptotic constraints (^) — (pD|). It follows that 
the appropriate sets of free functions, restricted to satisfy the asymptotic 
constraints, parameterize the corresponding families of AVTD polarized and 
half-polarized solutions. 

The work done in Sections 2-4 presumes a fixed choice of harmonic time 
foliation. In Section 5, we discuss the relation between alternative harmonic 
time foliations, and we show that one sees AVTD behavior using any one of 
them. The asymptotic relation between alternative harmonic time foliations 
explains why this makes sense, and hints at criteria for predicting whether 
one sees AVTD behavior using other choices of time foliation and observers. 



We show in Section 6, using the SL{2,R) Geroch transformation |19 
that there are U{1) Symmetric solutions, beyond those discussed here, in 
which AVTD behavior is seen. We speculate on further classes of solutions 
which exhibit AVTD behavior, and make concluding remarks, in Section 7. 
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2 ^(1) Symmetric Vacuum Spacetimes 



If we choose 27r-periodic coordinates ( ^) on T^, with the 

Kilhng field expressed as then the generic U{1) Symmetric spacetime 
metric on x may be written in the form 

g = [-N^e-^^dT^ + e-2-e^e«fe dx'^dx^] + e^'*'{dx'^ + padx^f , (1) 

with the functions 0, N , and A, and the form field (3a-, and the unit determi- 
nant symmetric tensor field Cab-, all functions of x"" and the time coordinate 
r (independent of x^). Here the latin indices a, h take the values 1, 2; note 
that we shall find it convenient to write = {u,v). 

The metric form (P does involve some restrictions on the gauge freedom, 
in that we have aligned the Killing field with x^ for all time, and we have 
set the corresponding shift vector to zero. We now further restrict the gauge 
freedom by setting the lapse function equal to e^. This choice of the 
lapse does not restrict the initial choice of Cauchy slice; however, once a first 
slice is chosen, all others are determined. We note that it follows from the 
condition N = that the time function r satisfies the equation 

□r = (2) 

where □ is the wave operator corresponding to the metric g.^ Hence this 
choice of slicing is sometimes called a "harmonic time."0 

Before writing down the full system of Einstein's equations for the U{1) 
symmetric spacetimes in terms of the metric (|l]), we wish to make a small 
restriction which considerably simplifies the analysis. To state this condition, 
we first write out the projection of the Einstein supermomentum constraint 
along the Killing field it takes the form Q 

r,a=0 (3) 

■^One has Dr — for □ corresponding to the 2 + 1 metric g ~ — iV^e^'^'^dr^ + 
e~'^'^ e^eabdx'^dx^ as weh as for □ corresponding to the 3 + 1 metric in equation (Q). 

"^This terminology is a consequence of the somewhat loose practice of using the word 
"harmonic" to refer to solutions of the equation Dr = as well as to the equation V'^r — 0. 

^The form of this constraint is reminiscent of the Gauss law constraint for electromag- 
netism. This reflects the fact that one may view 3 + 1 vacuum gravity with U{1) symmetry 
as a Kaluza-Klein theory for 2 + 1 Einstein-Maxwell fields, plus a Jordan scalar field. 
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where /° is the momentum conjugate to Pa- The general solution to can 
be written as 

r = e"^w;,,+/i'^ (4) 

where w is a scalar function, e"'' := ^ ) ' dual to a harmonic 

one-form on T^. For the x R spacetimes of interest here^ one finds (using 
the rest of Einstein's equations), that (/i") = ( ^2 j > fo^ constant 

in time as well as in space. 

So far, no restrictive assumptions have been made. Now, however, we 
set both and to zero. Correspondingly, we write j3a in terms of an 
integral involving one free function, which we call r (see for details), 
and it follows that the canonical pair (/?„, /°) is replaced by the pair (tu,r). 
These restrictions are consistent with Einstein's equations; they are some- 
what analogous to the vanishing of the two "twist constants" which defines 
the Gowdy spacetimes as a subfamily of the Symmetric spacetimes. One 
might wish elsewhere to consider what happens if or is nonzero, but 
we shall not pursue the issue here. 

If we now choose the following convenient (and non-restrictive) form for 
the unit determinant tensor Cab, 

_ 1 / + e-2^(l + xf e^^ + e'^'ix^ - I) \ , . 

2\e^-^ + e-^'{x^ -I) e^' + e-^'{l-xf J 

(where z and x are functions of u, v and r), then the Einstein evolution equa- 
tions for these spacetimes can be expressed as a Hamiltonian system for the 
fields {(f),A,w,x,z) and their respective conjugate momenta {p,PAr,Px,Pz}- 
The Hamiltonian for this system is 

H = h dudv (6) 
Jr2 



with 

1 2 

2" ' 8" ' 2^ ■ 2 



h=oPl + n^^'Pl + oP^ + o^^'^r^ - -pI + 2p 



^We note that if one studies U{1) Symmetric spacetimes on (E^ x S^) x R, for more 
general surfaces S^, the same equation appears, and one has the same sort of solution 
(^. However, for = S*^, (/i") necessarily vanishes, while for higher genus S^, {h"") can 
be more complicated. 
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+e-2- {(eV") (eV") ,aA,6+e^ [(e'^^) (e'^^) 



where e"'' is the matrix inverse of Cab (see equation (|)). The evolution equa- 
tions for {(p, A, w, X, z] p, Pa, r, p^, Pz} are obtained via Hamilton's equations, 
using this function H. 

There are also constraint equations which this data must satisfy. They 
take the form 

= Ho 

= h-2pA (8) 

o = nu 

= PzZ,u +PxX,u +PaA,« -Pa,u +P(f),u +rw,u 
+ ^{[e^'-{l + xy] p^ - {1 + x)pz} 

A{[e^^+{l-x')]p^-xpz},u (9) 

o = n^ 

= pzZ,v +PxX,v +Pa-^,v -PA,v +p(p,v +rw,^ 
~{[e''~{l-xf] Px + {l-x)pz},u 

+ ^{[e^' + {l-x^)]p,-xpz},y (10) 

The first of these is the superHamiltonian constraint. The others are super- 
momentum constraints. Recall that one of the supermomentum constraints 
-equation (^-has been solved, and is therefore of no further interest. For 
later purposes, it is important to note that Tio is written as a double density, 
while Tiu and Ti^ are expressed as single densities. This fact is irrelevant 
in seeking solutions to the constraints, but it is important in verifying the 
preservation of the constraints under evolution (see Section 4). 

In an appropriate sense, the Hamiltonian system (|^)-(0) with constraints 
(|^)-([T0|) (and with the remaining gauge freedom described above) has two 
gravitational degrees of freedom. If one can set one of the functions or w, 
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together with the corresponding conjugate momentum p or r, to zero on some 
Cauchy surface, and if the evolution equations preserve the vanishing of these 
two conjugate variables, then one has reduced the system to one gravitational 
degree of freedom, and the resulting family of solutions is said to be "polar- 
ized." Inspecting the Hamiltonian one finds that the only conjugate 
pair for which this can be done is (w, r). Thus, the polarized U{1) Symmet- 
ric spacetimes are those with the canonical variables {0, A, x, z; p, pa, Px, Pz}, 
with the Hamiltonian given by (P)-(0) with w = and r = 0, and with the 
constraints given by (|8|)- ([10|) with w = and r = 0. We note that for these 
polarized solutions, the spacetime metric takes the form (|1|), with j3a = 0. 

As noted in the introduction, we find AVTD behavior in a family of 
"half-polarized" U{1) Symmetric solutions, as well as in the family of po- 
larized U{1) Symmetric solutions as just described. For these half-polarized 
spacetimes, rather than setting a conjugate pair of variables to zero, one ties 
the conjugate pair together, so that effectively the pair involves one free func- 
tion instead of two. It is not known whether this can be done consistently in 
terms of initial data on a Cauchy surface. It can, however, be done in terms 
of asymptotic data, as we show in the next section. 

3 Fuchsian Study of the Evolution Equations 

While AVTD behavior can be described in terms of asymptotic approach to 
Kasner evolution of the metric fields at each spatial point (see Introduction), 
for the purposes of proof, the following formulation is more useful. 

Definition (AVTD Behavior) 

A cosmological solution {M,g) exhibits AVTD behavior if there exists a global 
spacelike foliation of{M,g), and there exists a solution {M,g) of the Ein- 
stein VTD equations (which are obtained by dropping all spatial derivatives 
in the evolution equations written with respect to St, and dropping all spatial 
derivatives in the super Hamiltonian constraint) such that if tg G [— oo, +oo] 
labels the singularity in {M,g), ther^ limt-,t^ \g(t) — g(t)\ = 0. 

^Although this definition in principle makes sense for any choice of norm | |, for our 
purposes here, we presume that | | is the absolute value norm for each of the components 
of g, evaluated independently for each spatial point, so the convergence is pointwise, L°°. 
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Thus, to show that the spacetimes in a given family of solutions have 
AVTD behavior, one first needs to find an appropriate foliation, with the 
singularity occurring at a well-defined value of the time parameter. While it is 
often useful in Fuchsian analyses to choose coordinates so that the singularity 
occurs at t = 0, for the U{1) Symmetric spacetimes, it is convenient to work 
with the time coordinate r (introduced in Section 2), in terms of which the 
singularity occurs at r ^ oo. Note that in terms of t = e~'^, the singularity 
occurs at t = 0; however the use of r instead of t simplifies the analysis here. 

We next identify the VTD equations for the U{1) Symmetric spacetimes. 
The VTD evolution equations^, written with respect to r and the coordinates 
chosen in Section 2, can be obtained by varying the Hamiltonian 



H 



rp2 



dudv . (11) 



As discussed in Section 4, the general solution to the VTD evolution equa- 
tions can be written out explicitly. In fact, for the purpose of showing that 
the polarized and the half-polarized U{1) Symmetric solutions have AVTD 
behavior, we only need a large r expression for the VTD solution. For the 
polarized case (r = 0,w = 0), we have the following: 

Large r VTD Solution (Polarized) 

^ o 

(t){u,v,T) =(j) {u,v) - v^{u,v)t (12) 
p{u, V, t) = -4v^{u, v) (13) 

o 

A{u,v,t) =A {u,v) + 2t — va{u,v)t (14) 
Pa{u,v,t) = va{u,v) (15) 

(16) 



®In this section, we ignore the constraint equations, both for the VTD system and for 
the full Einstein system. 



10 



(17) 



z{u, V, t) =z (m, v) — Vz{u, v)t (18) 

pz{u, V, t) = -4:Vz{u, v) (19) 

Here and are strictly positive f\, and the notation (p, p, etc, is used to 
indicate that these are VTD solutions, not Einstein solutions. The eight 

° o o o 

functions {0, f^. A, f a, 2;, f^, z, f^j are all - until we consider the VTD con- 
straints, in Section 4 - treated as free functions on T^; they parameterize the 
set of solutions of the large r VTD equations corresponding to H in equation 

m 

We now write the polarized U{1) Symmetric fields each as the sum of 
the appropriate term from the large r VTD solution, plus a remainder term. 
Specifically, for eO, we have 

Expansion Ansatz (Polarized) 

= + (50 (20) 

p=p + e~'^5p (21) 

A = A + (5A (22) 

Pa = Pa + e~'^(5pA (23) 

X = x + 6x (24) 

Px =Px + e"'^5px (25) 

z = z + 5z (26) 



^This parametrization of the large t VTD solution requires this sign condition. 
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Pz=Pz + e ^'^Sp, 



(27) 



The aim now is to substitute this ansatz into the Einstein evolution equa- 
tions, and show that the set of remainder functions {50, 6p, - ■ ■ , Sp^}, collec- 
tively, satisfies an evolution equation in Fuchsian form. That is, if one uses 
to denote a vector whose components consist of {6(j), 6p, ■ ■ ■ , 6pz} and cer- 
tain of their spatial derivatives, then one wants to show that satisfies the 
Fuchsian equation^ |2l[ [g 

(9,^-A^ = e-'^"F(r,M,t;,^,V^) (28) 

where A is a matrix which is independent of r, and for which 
(T'^[:= exp{Aina)] is uniformly bounded for < a < 1; where /i is a strictly 
positive constant; and where F is continuous in r, is analytic in u,v and 
and satisfies an estimate of the form 

\f{t,u,v,^,v^) - F{T,u,v,e,ve)\ < c[|^-e| + |v^- ve|] , (29) 

for some constant C, provided that 6, V6, \E', and V\l' are bounded. ( Here 
"II" refers to the pointwise absolute value norm, summed over the relevant 
components of F,"^, etc.) 

For our present purposes, one readily verifies that it is sufficient to show 
that each of the functions 6ilj,6p, etc (generically labeled "g") satisfies an 
evolution equation of the form 

drq-i^q = J2 e"''' Vfc (r, M, ^, V^) (30) 

k 

where the sum is finite; where z/(-u, v) ^ 0; where the fj,k{u, v) are r-indepenl- 
dent functions, each bounded from below by a strictly positive constant; and 
where the functions fk are analytic in u,v, r, and V\l/, and are bounded 
by a polynomial in r. 

We shall now show explicitly that the evolution equations for 6(f) and 
6p take this form, and argue that the evolution equations for the other six 
functions SA, etc, do as well. 

^°ln terms of i = e^^ , equation (^8|) takes the more familiar Fuchsian form 
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We start with the evolution equation for 50. Based on the Einstein Hamil- 
tonian equation 



1 

4^ 



(31) 



which we derive from the Hamiltonian (|]) and (0), we obtain 



(32) 



where we recall that e is the strictly positive parameter appearing in the 
expansion ansatz (p0| ) -(l27| ) . Clearly this is of the right form, with z/ = 0, 
jii = e, and /i = 6p. 

The equation for Sp is not so simple. Based on the Einstein Hamiltonian 
equation 



drP = 4e-'^ (eVV,a) 



we derive 

Ot-Sp — e6p = 4 



(33) 



(34) 



ith e the inverse of Cab from (g); i.e., 



ab 



^2z ^ ^-iz^l _ ^ 



-e2- _ - 1) e2^ + e-2^(l + x)2 



and with x and z to be expanded as in equations (^ 
Let us now write e"^ in the form 



(35) 

and f|T2D-f|T9D. 



-2z ■ 



[X 



1+X 



1 -1 
-1 1 



(36) 



The right hand side of (^^ splits, and we have 



drP — e6p = 4 < e 



^A+€T-vjyT+5A -2z+2vzT-25z ab 



+ 4 gA+eT-s;A'r+<5Ag+2z-2i;zT+2(52;ga6 



a -V^^aT + S(j)„ 



(37) 
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Calculating out the first term on the right hand side of equation (P7|), we 
obtain 

RHSj = e^^-A+2v.)r^J.-2l^SK~2S. I^ab^^ ^^^^ ^ ^^^^ ^ 

+ef (^0,ab -V<t>,abr + 6(p,ab 

+ A,b -2 -K,^ -2t;^,Jr + 5A,b -25z,b 

ef (^^,a-v^,ar + S^,a^Y (38) 

This expression contains second derivatives of some of the fields included in 
and so it cannot satisfy the conditions demanded for fk from equation 
(pop as listed above. However, by introducing a new variable 



0a = (P,a (39) 

together with its corresponding expansion with remainder term, 

o 

0a =0,a -V<t>,a'^ + ^a, (40) 

and by including 5(f) a in an expanded version of ^ , one readily transforms 
( p7| ) into a function involving only first derivatives of the expanded \E', as 
required. 



Now, inspecting the terms in (|38|) , including those in the matrix , we 



see that, so long as the functions v\ and satisfy the condition 

VK >2v, + e (41) 

for e > as introduced above (equations (pOD-(p7|), then RHSj takes the 
form required for the right hand side of (^), with /i/ := f a — — e. 
Calculating the second term on the right hand side of (pTf), we get a very 
similar expression, only now we have ixn := + — e and the condition 
for flu to be strictly positive is 

va > -2vz + e . (42) 

Assuming that Vz is strictly positive (which is a requirement for this expan- 
sion ansatz), (^) is the more rigorous condition. Thus we find that in this 
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first step of checking whether solutions of the ansatz form (PD[)-(P7[) are to 
be AVTD, it is sufficient that condition (|41|) hold. 

We proceed to consider the rest of the evolution equations for 6A, 6pA, etc 
(including the added first derivative variables). We get very similar results. 
For example, the 6A evolution equation is 

dJA = -e~'^6pA, (43) 

while the Sp\ equation - derived from the evolution equation for pa, 

- + 2eVV,a0,6 } (44) 



A,ab -VA,abr + 6 A, 



ab 



-can be written as 
dr6pA-e6pA = _e(^-''A+2''-)-e^-2Sg5A-25. jgafc 

+ef (Ji,a -VA^aT + SA,a ) (^A,b -VAfiT + 5A,f, 
+ 2ef (^,a -V^^aT + S(p,a ^ f^,^ -V^f^^bT + S(f),b 

+ similar terms} 

_g(e-.A-2.,)rgA+2Sg5A+25. j^^ab + . . . j . (45) 

We thus find that for this pair, too, so long as the inequality ( ^T]) holds, and 
so long as we define new field variables 

Aa = A,a , (46) 

•^a = -^la 1 (47) 

Za = Z,a , (48) 

with corresponding expansions 

Aa=A,a-VA,ar + 6Aa, (49) 
Xa =X,a +SXa, (50) 
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Za=Z,a-V^,ar + 5Za, (51) 

the evolution equations are of the proper (Fuchsian) form. The same inequal- 
ity condition (|iT|), together with the condition 

V, > e/4 (52) 

leads to the proper form for the evolution equations for the quantities Sz, Spz, 
Sx, and 5px, as well as for 6(f)a, S-^a, Sxa, and 6za- The evolution equations for 
these latter four quantities are derived /^from their definitions. For example, 
the evolution equation for S4>a is 

dr5(l)a = \e-'^da5p. (53) 
We have now verified the following: 



Proposition 1 If the polarized U{1) Symmetric gravitational variables 

Qpol ■= {(f),(f)a,P,-^,K,PA,X,Xa,Px,Z,Za,Pz} 



are expanded as in equations (\l^)-^T^), ^2di)-^li) , (U^), and then 
for any choice of the asymptotic data 



•Apoi ■= |0, A, X, z, v^, va, Vx, Vz 
which satisfies the conditions 

Vz > e/4 (54) 

and 

VA > 2vz + e, (55) 

(for e > 0), together with the expansion ansatz condition 

> 0, (56) 

the vacuum Einstein evolution equations take the form of a Fuchsian system 
for 

SQpoi ■= 50a, Sp, 5A, 6Aa, SpA, Sx, 5xa, Sp^, 5z, 5za, Spz} . 
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It then follows from the work of Kichenassamy and Rendall |^ that we have 



Corollary 1 For each analytic choice of the asymptotic data Apoi which sat- 
isfies conditions ^^-1^^, there is a unique analytic solution 6Qpoi to the 
vacuum Einstein evolution equations for sufficiently large r, with 5Qpoi ap- 
proaching zero as t approaches infinity. 

Hence we have an .4porparameterized family of asymptotically velocity 
term dominated solutions to the U{1) symmetric vacuum Einstein evolution 
equations. 

We have so far ignored the Einstein constraint equations. We shall con- 
sider them in Section 4. In the rest of this section, we show that there 
is another asymptotic expansion ansatz which allows for nonzero r and w 
and yet still shows AVTD behavior. It includes the polarized ansatz as a 
subcase; and it introduces one extra free function into the asymptotic data, 
rather than two. Hence we call these "half-polarized" solutions. 

The half-polarized solutions are defined by their asymptotic form, which 
appends the following expansions to those of the polarized solutions given in 

mvm) and mvm 



Expansion Ansatz ( Half-Polarized) 

r('U, t>, r) =r (m, f ) + e^^'^^r (57) 



(58) 

These expansions introduce one new function, f (m, f), into the set A1/2 
of asymptotic data, and include two time-dependent expansion functions 
5w{u, V, t) and 5r{u, v, r). 

To verify that Proposition 1 and its Corollary 1 extend to the half- 
polarized ansatz, we first expand out the evolution equations for 6w and 
6r. For 6w, it follows from the Einstein Hamiltonian evolution equation 

d^w = e^'^r (59) 



w{u, V, t) 



4 4>-v^T+S<f) 



+ 5w 
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that we have 



d^5w — Av^dw = e 



4:V, 



— 5w ] 5p + 5r 



(60) 



which is appropriate Fuchsian form, so long as f ^ > (a necessary condition 
for the expansion ansatz). For 6r, it follows from the Einstein Hamiltonian 
evolution equation 



-2t 



(61) 



that so long as condition (^) is satisfied, and so long as we introduce first 
derivative field variables for w, we have Fuchsian form as well. 

We still need to verify that the evolution equations for 50, 6p, and the 
other remaining functions in SQ1/2 retain Fuchsian form. The evolution equa- 
tions for 6(f), 6A, 6x, and 6z, as well as for their first derivative field variables 
50a, SXa, 6xa and 6za are unchanged by the nonvanishing of the r and w 
terms in the Hamiltonian (§)-(0), so they are in proper form. The evolu- 
tion equations for 6p, 6pA, 6px, and 6pz all pick up extra terms which need 
to be checked. For example, the equation for 6p, derived from the Einstein 
Hamiltonian evolution equation 



drP = 4e-2-(eVV,a) 



2e^V + 2e-2^e^e-^Vw;,aW;,b 



(62) 



takes the form 

drSp — e6p = 4 



_j_2gA+{':-'"A)T+(5Ag-40+4?)^T-4<50gafe ^ 



^i(l)-Av^T+AS(t> 



4:V, 



4:V, 



+ 6w 



+ 6w 



(63) 



The first term on the right hand side of equation (p3| ) is just the right 
hand side of equation (|34D , the evolution equation for 6p in the polarized 
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case; it has already been checked. The second term, a new one, clearly is in 
Fuchsian form so long as 



The third term can be written in the form 



(64) 



X < e 



ab 



+ 6w\ + 



+ 



4t;' 



X 



r 



+ + 



^,6 



_^^(e-VA-2v^-4:V4,)T^A+2°z+4^^5A+25z+4:5(P ^ 

X {cjj [same terms]^ [same terms]^} (65) 



So we find that this term-and consequently the half-polarized evolution 
equation for 5p-takes Fuchsian form so long as 



(66) 

which follows from conditions (|5^)-(|56D. Examining the extra terms in the 
evolution equations for 6pA, Spx, and 6pz, we reach the same conclusion. The 
full system of evolution equations (including the first derivative variables 
introduced above) takes Fuchsian form, and we have the following: 

Proposition 2 If the U{1) Symmetric gravitational variables 

Ql/2 ■= {(l),(pa,P,-^,K,PA,X,Xa,Px,Z,Za,Pz,W,Wa,r} 



are expanded as in equations ^^)-^%), then for any choice of the asymptotic 
data 



° ° o o o 

(p,A,x,z,v^,VA,Vx,Vz,r 



which satisfies conditions j^)-^^), the vacuum Einstein evolution equations 
take the form of a Fuchsian system for 



SGi/2 ■= {S4>, S4>a, Sp, SA, 6Aa, 6pA, Sx, 6xa, Sp^, Sz, 6za, 6pz, Sw, 6wa, Sr} 
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Corollary 2 For each analytic choice of the asymptotic data A1/2 which 
satisfies (|5^j-(|5^j, there is a unique solution 5^i/2 to the vacuum Einstein 
evolution equations for large t, with 5Qi/2 approaching zero as r approaches 
infinity. 

4 The Constraint Equations 

The results of section 3 show that for the polarized as well as for the half- 
polarized families of U{1) Symmetric spacetimes, there is a function space of 
asymptotic data such that for any choice of data in that function space there 
is a unique corresponding solution of the U{1) Symmetric Einstein evolution 
equations which is AVTD and asymptotically approaches that choice of data. 
These results ignore the Einstein constraint equations; hence the spacetimes 
they describe are generally not solutions of the Einstein vacuum field equa- 
tions. We now wish to show that for both of these families, there is a set of 
constraint equations on the asymptotic data such that if a set of asymptotic 
data satisfies these asymptotic constraint equations, then the corresponding 
solution of the evolution equations satisfies the full constraint equations (and 
therefore satisfies the full Einstein vacuum field equations). In addition, we 
wish to show that, for each of the two families, the asymptotic constraint 
equations admit a set of solutions of the appropriate number of parameters. 
We may then conclude that we have a corresponding set of AVTD solutions 
of the Einstein vacuum field equations in each of the two cases. 

We note that, in addition to the constraint equations Ho = 0, Hu = 0, 
and Tiy = 0, which we treat here, there are new constraints on the initial 
data of the form 0^ = da(f) which are necessitated by the introduction of the 
new variables (pa, ^a, Xa, and Za , as discussed above. These readily translate 
into constraints on the asymptotic data, and are preserved by the evolution 
equations, so we shall not treat them here any further. 

We begin the analysis of the constraints TCo = 0, Hu = 0, and Hv = 
by deriving the corresponding asymptotic constraint equations. We do 
this for each family by substituting the appropriate expansion ansatz for 
the fields into the constraint equations, and then letting r ^ 00, with the 
consequent vanishing of all 5(field) terms and all terms consisting of decaying 
exponentials times fields bounded by polynomials in r. More specifically, 
starting with the polarized spacetimes, we substitute expressions (|20|)-(p7[) 
into (H), (H), and ([T0|) . Then, setting 5x, 5z, 50, SA and their derivatives 
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to zero, and setting e ^^"^ and like terms to zero as well, we obtain the 
foUowin^^. 

Asymptotic Constraint Equations 



'Ho 
2 



(67) 



-AVz X,u +Va A,„ -4^70 (f),u -f A,n 



1+ X] V:r + 4:[1+ X] vA ,v 



+ (1— X ] Vx + 4: X V;> 



(68) 



= Ti:^, 

o o o ° 

:= -4f2 z,^ +Vx x,^ +f A A,„ -4t;<^ -ua . 
+ ^|(l-^)'^^x + 4(l-x) t;, 

+ ^{(l-^')t^x + 4xi;,}„ (69) 

These are constraints on the choice of the asymptotic data Apoi for polarized 
f/(l) Symmetric spacetimes. A similar procedure leads to constraints on 
Ai/2, the family of half-polarized U{1) Symmetric spacetimes. One finds 
that the asymptotic constraints here are exactly (^) - (|69[) ; the terms ^e^'^r^ 
and ^e^e~'^'^e"''^w,a w,b in H and the terms rw,u in 7Y„ and rw,^ in have 
no asymptotic effect if one expands r and w as in (|57D and (0). 

We now want to argue that if the asymptotic data Apoi, or A1/2 satisfy 
the asymptotic constraints (|6^)-(|69|) and if the corresponding fields Qpoi or 
Qi/2 satisfy the U{1) Symmetric evolution equations generated by (^, then 

00 r- — . 

^^In calculating Jiu and Hv, we have used to cancel certain terms and thereby 
simplify the expressions and ( |69| ) 




21 



Qpol, and Q1/2 satisfy the initial value constraints (Q)-(|TU|) on any Cauchy 
surface sufficiently close to the singularity (i.e., for sufficiently large r). We 
argue this as follows (focusing first on the polarized case): If we substitute 
the expansion ansatz (pOD-(p7D for the polarized fields into the expressions 
([^)-([T0|) for the constraints, then we can write 

Ho =m +e-^^5Ho (70) 

rtu =nu +5nu (ri) 
K =nv +sn, (72) 

where {Ho,'Hu,'Hv} are the asymptotic constraint expressions (p7D-(|5UD, fi 
is a strictly positive constant, and {6H0, SHu, STCy} are the remainder terms 
(defined by (|70| ) - (|72|) ). Our calculations leading to (^)-(iD show that, re- 
gardless of whether the asymptotic data Apoi satisfy the constraints (|67D-(|69D, 
if the fields Qpoi evolve via the U{1) Symmetric evolution equations, then the 
remainder terms approach zero as r — 00. We want to show more. We will 
show that in fact, so long as the asymptotic constraints (|67|)- (|69|) are satis- 
fied, Hq, TCu, and Hv are identically zero for sufficiently large r. We do this 
by using the U{1) Symmetric evolution equations to show that {Hq, Hu, 'Hv} 
are identically zero for sufficiently large r. To do that, we use the U(l) Sym- 
metric evolution equations to show that 57^0; ^'Hu-, STi^ satisfy a Fuchsian 
system, and then note that this Fuchsian system admits 6Ti.o = 0, 6Hu = 0, 
STiv = as a solution. Thus, since the solution to the Fuchsian system (for a 
given set of asymptotic data) is unique near t 00 , and since {0, 0, 0} is a 
solution, it follows that we have SHq = 0, 6Hu = 0, STCy = for sufficiently 
large r. To complete the argument, we note that if the asymptotic data 
Apoi satisfy the asymptotic constraints and if {SHq, SHu, SHv} = {0, 0, 0} for 
large r, then it follows ihom (^-(^ that Ho = 0, Hu = 0, = for 
large r. The constraints are satisfied. 

The key to this argument is the derivation of the Fuchsian system for 
{STCo,6Hu,SHv}- To obtain this system, we start with the evolution equa- 
tions for {'Ho,'Hu,'Hv}, which can be calculated from the divergence Bianchi 
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identity:[] 

drHo = -e^-2-[A,,e'^''(7^, - m^) + (e'^'(K - /Sa^s)),;,] (73) 

drUa = -daHo (74) 

where a, b run over u and f . If we now substitute the expansions ([70|)-(|7^ 

o o o 

into these evolution equations, and further set Ho^Hu^Hv and equal to 
zero, then we derive 

dJHo - fiSHo = (e^-2^+'^^) [A,h e'^'SHa + {e'^'SHa) ,b ] (75) 

dJHa = ~e-^^da mo) (76) 

These equations clearly constitute a Fuchsian system for {SI-Cq, STCu, SHy} so 
long as /i > and f a > 2vz + /i, which are the familiar inequality conditions 
on the asymptotic data. Further, we note that {6T-Co,6Hu,SHv} = {0,0,0} 
is a solution of the system ([f5|)- ([76|) . Finally, we note that (|75D -(|76D holds 
not just for the polarized U{1) Symmetric spacetimes, but in fact for the 
half-polarized U{1) Symmetric spacetimes as well. Thus we have proven the 
following: 

Proposition 3 // the asymptotic data for either a polarized or a half-po- 
larized U(l) Symmetric spacetime is analytic and satisfies the asymptotic 
constraint equations j^^-j^), then the corresponding spacetime satisfies the 
full set of Einstein constraint equations-and hence the full set of Einstein 
vacuum field equations-for sufficiently large t. 

While Proposition 3 tells us that any solution of the asymptotic constraint 
equations (p7D-(|5UD leads to a solution^ of the vacuum field equations, it 
does not tell us anything about finding and parameterizing solutions of the 
asymptotic constraints. We could work directly with (^)- (l69|) ; however, 

^^Note that while Ho is a weight two scalar density, TCu, Hv and TCs are all weight one 
scalar densities. Note also that while we have already eliminated Ti,3 = (see equation 
(H)), we leave Tis in equations (|7^)-(|7^) for the sake of generality. 

^■^In fact, since ro and wq do not appear in the asymptotic constraints, every solution of 
(|67|)-(|69|) leads to a large set of half-polarized solutions, along with one polarized solution. 
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since they are specified on a possibly degenerate manifold, we choose an 
alternative approach. 

The alternative approach starts with the recognition that (|67D-(|69D are 
the asymptotic constraint equations for the U{1) Symmetric VTD equations 
as well as for the U{1) Symmetric vacuum Einstein equations. To see this, 
we substitute the large r VTD solutions-(|12])-(0) for the polarized case, and 
([T^)-(p^) together with r =r and w = for the half-polarized case-into the 
VTD constraint equations: 



= Ho 

= ip2 , I 4. 2 



-p H r Pf, 

8^ 2 2^^ 



(77) 



Hu (see equation (9)) 



(78) 



Hv (see equation (10)) 



(79) 



We obtain, in both cases, equations (p^)-(|69|). 

We next note that we have global existence in time r for solutions of 
the U{1) Symmetric VTD evolution equations. This property holds even for 
the slightly generalized form of these equations, which we obtain ^from the 
Hamiltonian 



H 



rp2 



dudv. 



(80) 



where a is a freely specifiable "lapse" function. Global existence is an imme- 
diate consequence of the following expression for the general solution to these 
evolution equations (clearly well behaved for all values of re(— 00, +cxd)) ||^ 



(ar - r^) - -in 



m (l + e" 



■4:Vj,i ar-T, 



(81) 



p = -Ava 



I _ g-4t^.^.(,"T-T-0J 



l + e" 



(82) 
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w 



1 + e"'"^ 



X 



Z = —V, 



(83) 

r = -AC^v^ (84) 

= -4C.t^. (86) 
(«r - r,) - ifn (l + e-^^^^^^"-))] (87) 



= -4-^ ( i^,-4..(..-..) ) (88) 



A = Ao + 2r - avA_r (89) 

PA = VK (90) 

Here r^, > 0, ^</), t^, "^^ > 0, Cz, ^t/-, Aq and are all free functions of 
u and f ; they are "constants of integration" for the set of ordinary differen- 
tial equations which comprise the U{1) Symmetric VTD evolution equation 
system. Note that (|8T|)- (|90|) is the general solution for the full U{1) symmet- 



ric VTD evolution system; no polarization or half-polarization assumption 
has been made. Note also that the positivity conditions on v^f, and Vz are 
necessary for this parametrization of the solutions of the VTD equations. 

With global existence of the VTD solutions established, we may consider 
solutions of the VTD constraints ([r^)- (|79|) at finite times r, at which there is 
no difficulty with degeneracy of the spatial geometry. To relate solutions of 
([77|)- ([79|) to solutions of the asymptotic VTD constraints (|67D-(|6^) (the same 



as the asymptotic Einstein constraints), it is useful to establish the following 

Proposition 4 The VTD constraint functions TIq, Ihiu, andiiy, when evalu- 
ated for solutions (|5ip-(j£^ of the VTD evolution equations, are independent 
of T, so long as for some Tq, Hq^Tq) = 
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Proof: The time independence of the quantity Tio may be estabhshed di- 
rectly, by substituting expressions (pTD-(pO|) for the fields {0(t), p{t), w{t), 
r{T), x{t), Px{t), z{t), Pz{t), A{t),pa{t)}, into the expression (^) for Hq. 
One obtains, for all r (and all a), 2v1 + 2t>^ — |t>^, which is independent of 

r (and equal to 'Hq). 

Rather than proceeding to substitute (|HT|)-(P(]|) into (^) and (|75D, we may 
establish the r independence of 7t!„ and H.^ using the following argument: 
Since Idu = 'Hu and 71^ = H^, we find that Hu and both generate spatial 
diff'eomorphisms, tangent to r = constant surfaces. With Hq constant and 
presumed zero at time tq, we have Hq = for all time. Hence for any vector 
field M"'[x) on T^, we have 

[ M''{x)na{x)dx, [ a{y)no{y)dy] = (91) 

where { , } indicates the Poisson bracket. 

The generalized VTD Hamiltonian H, from equation (|5(]|), generates VTD 
time evolution for 'Hu and H,^ (as well as for any other function of the VTD 
fields). As argued above, Hu and H.^ commute with / aHo. Calculating the 
Poisson bracket of H.u and H.^ (times an arbitrary M", integrated over T^), 
with the remaining piece of the Hamiltonian H, we have 



M''{x)na{x)dx , 2 / pK{y)dy 

rp2 Jrp2 



{My^),aMx,2 / pK{y)dy 

rp2 



= -2 / {MyA),adx 

= (92) 

Thus iiu and Ti^, commute with the VTD Hamiltonian, from which it follows 
that H,u and H,v are VTD constant. 

Our remaining task now is to show that one can find a full set of solutions 
of the VTD constraint equations (|77|)-(|79D - five|3 free functions on in the 
polarized case; and six free functions as in the half-polarized case. To 
do this, it is useful to first show the following: For any fixed choice of the 
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Three of these functions are related to coordinate choice, and thus are not physicaL 
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asymptotic data Apoi (or A1/2) there is a choice of the lapse function a for 
which the VTD fields at r = 1 have constant mean curvature (CMC) in a 2 
+ 1 dimensional sense; i.e., based on the 2 + 1 Lorentz signature metric 

7 = -a^e^^-^^dr^ + e-^^+^eabdx'^dx' . (93) 

(compare with the 3+1 metric (1).) The spatial volume element for 7 from 
(H) is 

V = e^~'" (94) 

Using the VTD evolution equations and solutions to calculate the mean ex- 
trinsic curvature for a r = constant surface in this 2 + 1 dimensional space- 
time, we obtain 



trK 



-1 g.(V) 

(lapse) 

gA~2r 
^^A 



gAo-OTDA ' 

If we now set r = 1 and trK = C > in (^), we can solve for a: 



(95) 



a = —(Aq{x) - £n{vA{x)/C)) (96) 

Noting that our main results concerning the evolution behavior of the fields 
(Corollaries 1,2, and 3) hold only if va is positive definite, we see that indeed, 
for any fixed data Apoi or Ar, we can choose a (as in (0)) so that the r = 1 
hypersurface has constant mean curvature. We ensure the positivity of a by 
a sufficiently large choice of the constant C. 

Now that we know that any VTD solution admits a hypersurface with 
constant mean curvature in the 2+1 sense described above, we can use this 
CMC condition to help in the analysis of the VTD constraint equations, and 
show that we have a full set of solutions. We carry out this analysis using 
the conformal method, adapted to U{1) Symmetric data in 2+1 form, as in 
pO| (where it is applied to the Einstein constraints). Specifically, on T^, the 
manifold transverse to the U(l) orbits, we choose for our conformal data (1) 
a Riemannian metric Xab, (2) a symmetric tensor a'^'^ which is divergence- free 
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and trace-free with respect to \ab, (3) a constant r representing the mean 
extrinsic curvature of the geometry on T^, and (4) a pair of functions (j) 
and P(j) which are treated as matter fields, but actually correspond to the 
geometry in the direction tangent to the orbits of the U{1) isometry group. 
We then attempt to solve a system of three equations for a vector field W 
(which generates the longitudinal part of the extrinsic curvature on T^) and 
a positive definite conformal factor ip. If indeed for the given set of conformal 
data {X,a,T,(j),Pip}, solutions W and exist, then we may construct from 
W and ijj and the conformal data a set of initial data {7, tt, 0,p<^}, which 
satisfies the constraints. 

The equations to be solved for W and are generally coupled, and it is 
generally not easy to determine if solutions exist. For CMC conformal data, 
however, the equations decouple; and in the case of the Einstein constraints, 
it has been determined that solutions to the equations exist for essentially 
all choices of the conformal data. [pO| 



For the VTD constraints, of interest here, this decoupling still holds. 
While there has not been any systematic study of the equations for W and 
ip resulting from the VTD constraints, such a study is easily carried out. 
Indeed, the supermomentum constraints are the same for the Einstein and 
the VTD cases; in both cases, we have a linear elliptic system to be solved 
for W, and in both solution exists for all choices of conformal data 

satisfying the integrability condition 

pY^Va4>) = (97) 



where Y is any conformal Killing field on T^. The superHamiltonian con- 
straint, to be solved for ip, is different for the Einstein and the VTD cases. 
In the former case, it is a nonlinear elliptic equation for ip, while in the latter 
case, it is an algebraic equation for ip. In both cases, one verifies that for 
nonzero constant mean curvature and for generic choices of the conformal 
data, a unique solution exists. We note that if the conformal data is cho- 
sen to be analytic, it follows from the ellipticity of the constraint equations 
that the corresponding initial data is analytic. It further follows from the 
VTD evolution equations and their general solution that the corresponding 
asymptotic data is analytic. 

The analysis which we have just sketched shows that we have a full set 
of analytic solutions of the VTD constraint equations, and hence have a 
full set of analytic solutions of the asymptotic Einstein constraints. If such 
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solutions are to correspond to AVTD solutions of the Einstein equations, 
we must verify that the asymptotic data satisfies the inequalities Vz > e/4, 
V(j) > e/4, and v\ > 2vz + e, which are sufficient for the evolution equations 
to be Fuschsian. 

General solutions of the asymptotic constraints do not satisfy these in- 
equalities. We claim, however, that there are open sets of choices of the 
conformal data which do guarantee that these inequalities hold. To argue 
this, we first note that as a consequence of the asymptotic constraint equa- 
tion (|67[) va^ = 4:{vz'^ + v^'^), of the compactness of T^, and of our freedom to 
choose e to be any positive constant, all of the necessary inequalities (|5^)-(|56[) 
hold (for some e > 0) so long as and Vz are both positive definite. 

Now, in finding data which satisfy the VTD constraints at finite time 
(r = 1), we have no apparent direct control over the values of the data 
for arbitrarily large r (ie, the values of the asymptotic data). However, 
examining the explicit form of the global-in-time VTD solutions (|8Tl - pOD , we 
find that with r = 1 and = = r^, we have 



p = — 4t>^tanh(2aW(/,) (98) 

and 

Pz = -4:Vz tanh(2af ;2) (99) 

It follows from (|98| ) - (|99|) that for any set of data {(f),A,r,x,z,p,p\,u!,Px,Pz} 
at r = 1 which has p < and pz < for all u, v in T^, there exists a unique 
choice of positive definite functions and Vz such that (il)-(iD hold. (The 
same is true for nonzero, but sufficiently small, values of and r^.) Thus 
each solution of the VTD constraints at r = 1 which has p < and pz < 
has corresponding asymptotic data which satisfies the inequalities (|5^-(|5BD. 

The conformal procedure for solving the VTD constraints, sketched above, 
allows us to freely choose the function p at r = 1, and hence its sign; however, 
since the quantity pz is part of the extrinsic curvature of the two dimensional 
geometry in our 2 + 1 treatment of the constraints, we cannot directly choose 
Pz- Yet we can argue as follows that we can control the sign of Pz'- As just 
noted, Pz is related to the extrinsic curvature of the geometry on T^. More 
explicitly, letting vr"^ denote the momentum representation of this extrinsic 
curvature, we have 

= [7r""(e2" - e-2^(l + x^) + 'K'"'{e^' - e-^\l - xf) + 

2vr™(e2^ - e~2-(x2 - l))]e^-^^ (100) 
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Now in solving the momentum constraints via the conformal method, one 
determines the vector field W, which partially determines vr'^*. However, one 
can always add to this part of vr"'' an arbitrary divergence- free, trace- free, 
symmetric tensor density a"^. Since all metrics on are conformally related 
to a flat metric, one flnds that in appropriate coordinates, a"^ is divergence- 
free and trace-free if and only if is a spatially constant trace-free tensor 
density. It follows from this fact, and from equation (99), that through the 
choice of o""*, one can always guarantee that pz is negative deflnite. We then 
obtain, gued above, the positive deflniteness of Vz 

We conclude from this argument that, while it is not true that all sets of 
VTD data satisfying the VTD constraints at time r = 1 lead to asymptotic 
data satisfying the inequalities, if we impose certain open conditions on the 
choice of the conformal data, then we do obtain VTD data whose asymptotic 
data satisfies the inequalities. Noting the invertibility of functions such as 
f{v) = — 4tanh(2Q;u) for / < and av > 0, we readily verify that for both 
the polarized and the half-polarized U{1) Symmetric spacetimes, we have 
a full set of asymptotic data-two free functions for the polarized solutions 
and three free functions for the half polarized solutions (after quotienting 
out the diffeomorphism gauge freedom on T^)-which satisfy the asymptotic 
constraints as well as the inequalities which we have found to be sufficient to 
guarantee that the evolution equations are Fuchsian. Combining this result 
with those of Section 3 we have 

Theorem 1 There is a family of U(l) Symmetric solutions of the vacuum 
Einstein equations on x R which is AVTD with respect to a harmonic 
time foliation, is characterized by analytic asymptotic dataAi/2 satisfying the 
asymptotic constraints j^^-j^), and is parametized by three free functions 



on T^. 

Those members of this family of spacetimes which have asymptotic data 
with r= are AVTD solutions of the polarized U(l) Symmetric vacuum Ein- 
stein equations. This polarized subfamily is parametrized by two free functions 
on T^. 
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5 AVTD Behavior and the Choice of Observ- 



ers 

The definition of AVTD beliavior wliicli we use fiere depends on a clioice of 
observers; or equivalently, on a clioice of coordinates. One may therefore ask 
if alternative sets of observers will agree on the presence or absence of AVTD 
behavior in a given spacetime. This issue has not yet been carefuly addressed, 
since in most previous studies [|| |TT|, the family of spacetimes under 
consideration and the analytic techniques used have singled out a particular 
choice of time foliation and time threading. For the present study of U{1) 
Symmetric spacetimes, however, our choice of time is a bit more flexible, so 
we may begin to address the issue of the dependence of the verification of 
AVTD behavior on the choice of spacetime observers. 

As noted earlier, we use "harmonic time" in working with the U{1) Sym- 
metric spacetimes, which means that the time function r satisfies the wave 
equation Dr = 0. This condition does not fix the choice of time foliation; 
we may freely choose an inital Cauchy surface, along with an initial choice 
of lapse function (This amounts to Cauchy data for the function r). The 
condition Dr = then determines the rest of the time foliation. 

Fixing the time foliation does not necessarily fix the timelike observers. 
They are determined by the the choice of the spacetime "threading" ; ie, the 
choice of a congruence of timelike paths (corresponding to the worldlines 
of the observers). Our analysis is considerably simplified, however, if we 
require that the observer paths be everywhere orthogonal to the leaves of the 
spacetime foliation (so the shift vector field is everywhere zero). With this 
condition imposed, the choice of spacetime observers is fixed by the choice 
of the (harmonic) time function r. 

While it would be useful to compare the observations relevant to AVTD 
behavior that are seen by each set of observers corresponding to each choice of 
a harmonic time in a given U{1) Symmetric spacetime, we shall here pursue 
a more modest goal: Given a fixed polarized U{1) Symmetric spacetime 
and a fixed choice of harmonic time such that the corresponding surface- 
orthogonal observers see AVTD behavior, we shall show that there is a full 
(two free functions on T^) family of other harmonic time choices such that 
their corresponding surface-orthogonal observers see AVTD behavior as well. 

We call the parametrization of this family of time choices by two free functions " full" 
since the space of Cauchy data for solutions of Dt = consists of two functions. 
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Note that while we focus on the polarized case here to simplify the discussion, 
the same sort of results hold for half-polarized U{1) Symmetric solutions as 
well. 

We start by identifying the family of alternate harmonic time choices. 
To do this, we fix a polarized U{1) Symmetric spacetime (T^ x R,g) and a 
harmonic time choice r whose corresponding observers see AVTD behavior. 
We then write the harmonic time condition for a new time function T; in 
first order form, with ( := drT, we have 

drT = C (101) 



9.C = (e-'"eVT,J,,. (102) 

(Here we are choosing to work with the 2 + 1 rather than 3 + 1 version of the 
harmonic time coordinate). Now noting that 

f = a{x) + 6(x)r (103) 



C = b{x) 



(104) 



for arbitrary functions a{x) and b{x) > on solves the VTD version of 
equations ( |101| )-( p!0^ ) (which sets the right hand side of equation (|102|) to 



zero), we seek a family of solutions of ( |101| )-( piI2D of the form 

T = a{x) + b{x)T + 6T 



(105) 



C = b{x)+e-^^6C (106) 

for some e > 0. (The condition 6 > is required so that T uniformly 
approaches infinity as r approaches infinity). Plugging these forms into ( 101 )- 
( |102| ), we derive 

dJT = e~'^5C (107) 

and 

drSC - e5C = e(^-^-+2^^)V(r, w,t',T,Tv), (108) 
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where Ty denotes a set of new variables which we introduce so that only first 
derivatives appear, and where the function / is analytic in {u, v, r, T, Ty) and 
bounded by a polynomial in r. Thus, presuming the usual restrictions on va 
and Vz, we have a Fuchsian system for 6T and 6( such that T and ( of the 
form ( |105D and (|106|) solve (|101D -( [T0^ ). This gives us our family (parametized 
by the two functions a{u,v) and b{u,v) > 0) of choices of harmonic time for 
our fixed spacetime {T^ x R,g). 

Before addressing the question of observer-dependence of AVTD behavior 
for these different foliations, we wish to compare the paths of the observers 
orthogonal to the T = a + br + 6T foliation with those orthogonal to the r 
foliation. One way to do this is to start with the metric g{x^, r) expressed in 
terms of the r foliation and the observers orthogonal to the r foliation, then 
re-express g in terms of the T foliation but with the r-compatible observers 
retained, and finally determine the shift vector M{x,T), which measures the 
extent to which the r-compatible observers fail to be orthogonal to the T 
foliation. To carry out this calculation in practice, what we do is invert the 
transformation ( |105D , obtaining Q 

r = {-a{x)/b{x)) + {l/b{x))T-{l/b{x))6T{x,T) 

= c{x) + h{x)T + 5T{x,T), (109) 

and then substitute both r(x, T) and 



dr = {dkC + Tdkh + dk6T)dx^ + {h + ^^)dT (110) 

into the expression (^ for g{x,T), which for convenience we rewrite as 

g{x, t) = —n^i^x, T)dT'^ + •jijdx^dx^ . (HI) 

Here the indices i,j, k run from 1 to 3; and n{x, r) and 'Jijix, r) are certain 
combinations of e^, e^, etc (see (P along with the discussion just below it). 
We obtain 

g{x, T) = -n^h + dT^TYdT"^ - 2n^{h + dT5T){diC + Td,h + di5r)dTdx' 
+ - "-^(^ic + Tdih + di5T){djC + Tdjh + dj5T)]dx'dx^ 

(112) 

^^Our assumption that & > and our verification that 5T decays to zero for large r 
guarantees that there exists c(x) = {-^^), h{x) ~ > 0, and 5T{x^t) ~ with 



liniT^oo <5t = such that (|109| ) is indeed the x-parametrized inverse to (105) 
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where n{x, T) = n{x, c{x) + h{x)T+ST{x, T)), etc. If we now compare expres- 
sion ( |112| ) for g with the standard lapse-shift ( "ADM" ) form of a spacetime 
metric 

g = -N'^dt^ + hij{dx' + M'dt){dx^ + M^dt), (113) 
then we find that the shift vector M can be expressed as 

M' = -y^[n\h + dTST){djC + Tdjh + dj6T)] (114) 
where A*-' is the inverse of the tensor appearing as the coefficients of dx^dx^ 



m 



To see the behavior of the shift M* for large values of T, we need to write 
out n and X^^ in terms of the variables A, 0, x, z, and then substitute into the 
asymptotic expressions for these quantities. In doing so, we again need to 
replace r by its expression (|109|) in terms of T, noting that such quantities 
as 6 A which decay for large r do so as well for large T. We find that 

M' ^ (p^^-^K)hT ^ (polynomial inT). (115) 



Thus, recalling our usual restrictions (|54|)-(|56f) on t>A and f^, and also recalling 
that > 0, we see that the shift decays rapidly to zero for large T. It is 
especially telling that the shift decays to zero much more quickly than does 
the lapse as one approaches the singularity, since we calculate 

M'Mi/N^ ^ g(2^;,-,;A)/xT ^ (polynomial inT). (116) 

So we conclude that, for any values of the harmonic time parameter functions 
a{x) and b{x) > 0, the paths of the surface orthogonal observers correspond- 
ing to T = a -\- br + 6T become increasingly parallel to those of the surface 
orthogonal observers corresponding to r. 

We now discuss how to verify that the observers which are orthogonal to 
the T = a+br+ST foliation see AVTD behavior in the chosen spacetime. The 
key first step is to find coordinates y'^ which correspond to the T orthogonal 
observers. In practice, what we seek is a coordinate transformation 

(x,T)-^(x(y,T),r) (117) 

which removes the shift term; ie, for which g has no dydT term. A bit of 
calculation shows that if we choose the function x{y,T) so that it satisfies 

dTx'{y,T) = -M'{x{y,T),T), (118) 
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where M* is the shift vector field discussed above, then indeed the metric 
written in terms of the coordinate {y,T) has vanishing shift, and so indeed 
the y = constant observers are T surface compatible. 

We note three features of the coordinate y generated via ( |118|) . First, the 
ODE system ( |118|) together with a set of initial conditions Tq) = Ciy) 
constitute a well-posed system, with a unique local solution. Second, al- 
though the system (|118|) is generally nonlinear, it follows from the bounded- 
ness of M* that (|118|) (together with the chosen initial conditions) determines 
a unique solution T) for all T. Third, since M* decays to zero exponen- 
tially rapidly for large T, we may write solutions x\y,T) in the form 

x\y,T)=X\y)+6x\y,T) (119) 

where Sx^{y,T) decays to zero for large T. 

With the transformation between x and y determined, we next rewrite 
the metric g{x,t) in terms of the new coordinates. We are only interested in 
this for large T, so we may use expression ( |119|) , /^from which we derive 



dx' ^ X\j dy^ - M'dT. (120) 



Substituting this expression together with ( |119|) into formula ( |112| ), and 
also replacing n and 7,^ by the appropriate expressions in terms of Apoi = 
{(po, Ao, Xo, Zo, f<^, f A, Vx, Vz} and T, we obtain a coordinate representation of 
the metric of the form 

g{y,T) = gTT{y,T)dT' + gi,{y,T)dyW, (121) 

with gTT and gij as fairly complicated functions involving Apoi along with 
the functions c, h, and X\y). 

To show that the T surface orthogonal observers see AVTD behavior, it 
is sufficient to verify the following: 

Proposition 5 For every choice of the asymptotic data Apoi md for every 
choice of the harmonic time transformation functions {a, b > 0}, there exist 
functions {(^(^y),Ao{y),Xo{y),Zoiy),v^{y),VA{y),v^{y),Vz{y)}, and functions 
{S(j){y, T), 6A{y, T), 6x{y, T), 6z{y, T), 6p^{y, T), 6pA{y, T), 6p^{y, T), Sp^iy, T)} 
decaying to zero for large T, such that the metric coefficients gxr and gij in 
( \121\ ) can be written in the polarized U{1) Symmetric AVTD form described 
in sections 2 and 3; in particular, one has the asymptotic behaviors 
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^7.-7™n^^ (123) 

The proof of this proposition is a straightforward (somewhat tedious) 
consequence of carrying through the details of calculation ( |121|) . We omit it 
here, together with the explicit (not especially useful) expressions one obtains 

for {0(?/), Ao(i/),Xo(?/),5o(y),{;0(y),{;A(l/),^x(l/),^2(?/)} in terms of Apoi and 
{a,6 > 0}. 

Do we see AVTD behavior in these spacetimes using other foliations and 
sets of observers? This is not yet known. The foliations and observers we 
have discussed here exhibit two important features asymptotically: 1) All of 
the sets of surface compatible observers become parallel as one approaches 
the singularity. 2) In all of them, one sees AVTD behavior. It is not clear 
whether these features are related or not. It may be, for example, that in a 
spacetime which shows AVTD behavior with respect to one foliation and set 
of observers, AVTD behavior is seen by the surface compatible observers of 
any other foliation if and only if those observers become parallel asymptot- 
ically to the original ones. Or, it might be the case that these two features 
which coincide for our harmonic foliations and observers in these U{1) Sym- 
metric solutions are not closely related more generally. In either case, it may 
be that the surface compatible observers corresponding to most if not all 
spacelike foliations in a spacetime with AVTD behavior are asymptotically 
parallel, as one approaches the singularity. We hope to explore this issue in 
future work. 



6 Generating Further U(l) Symmetric Space- 
times with AVTD Behavior 

Numerical work |^ suggests that general U{1) Symmetric vacuum solutions 
on X i? show Mixmaster behavior near the singularity. We may then 
ask if the U{1) Symmetric solutions with AVTD behavior extends beyond 
the classes we have discussed here so far. In this section, we use Geroch's 
SL[2, R) method of generating new U{1) Symmetric solutions from old ones 
to show that this is the case. 

We recall |^ how the Geroch transformation works: If {0, A, w, x, z; 



P<t>,PA,r,Px-Pz} is a solution of the U{1) Symmetric vacuum Einstein equa- 
tions on (Ti^, R), and if ^ ^ is a (constant) matrix contained in SL{2, R) 
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( so det^ ^ ^ ^ = 1), then {(l),A,w,x,z;p^,pA,f,p^.pz} with 



e'' - ...... , ' , ...... , (124) 



c2(w2 + e4<^) + 2cdw + d? 



^ _ ac{w^ + e^) + {ad + bc)w + bd 

^ ~ c2(«;2 + e40) + 2cd«; + d2 ^ > 



^"^ c''{w^ + e^^) + 2cdw + d^ ^ ^ 



f = p(c^w + cd) + r[d^ + c2(w^ - e^^) + 2cdw] (127) 

is also a solution of the U{1) Symmetric vacuum Einstein equations on x i?, 
where is diffeomorphic to if and only if J^^g f = /^.g r. 

We wish to consider Geroch transformations which map solutions on E^ = 
to others also on T^. Since J^g r = if and only if = T^, we seek 
solutions and transformations such that J^,., r = both before and after the 
the transformation. Using the three constants of the motion 

! {2wr+p^), (128) 

JS3 



r, (129) 

and 

C= / [r{e^'t>-w^)-p^w], (130) 

Je3 

one finds that if {A,B,C} characterize the original solution, then 

A = {ad + bc)A + 2bdB - 2acC, (131) 

B^d^B- c^C + cdA, (132) 
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and 



C = c?C - \?B - abA (133) 

characterize the transformed solution. So we seek solutions and SL{2, R) 
matrices f ^, I such that B = and d^B — c^C + cdA = 0. 



c d 

For polarized solutions, r = = w and ^ 0, so B = 0, C = 0, and 
A ^ 0. Thus, to obtain B = 0, we need to choose ^ ^ ^ ^ with either 

c = OT d = (we can not have both zero, since ( ^ ] G SL{2,R). 



^ c d 

). For half-polarized solutions as well, the asymptotic decay of w and e^"^ 
require that C = 0, so again we obtain i? = if and only if either c = or 
d = 0. 

Applying an SL{2, R) Geroch transformation with c = to either a po- 
larized or half-polarized solution produces nothing really new. |^ The d = 
transformation, however, produces geometrically new AVTD solutions. Ap- 
plied to a polarized solution, such a transformation takes the form 

= {l/c^)e-^^ (134) 

w = a/c, (135) 

Pc^ = -P4> (136) 

and 

f = 0. (137) 

The new solutions are geometrically distinct from the original ones, since 
J e^'^dx^ is the diameter of the (evolving) three-geometry along the U{1) 
symmetry direction; and while this diameter decays to zero in the original 
solution, it blows up in the transformed solution. Note that one readily 

^^Such a transformation results in a (constant) rescaling of some of the quantities, and 
the addition of a constant to w. We recall that w does not appear in the metric, and it 
appears in the Hamiltonian and in the field equations only as a gradient term Vw. 
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verifies that so long as the original solutions are AVTD, the transformed one 
are as well. 

[ah 

For the half-polarized solutions, the transformation corresponding to 



takes the form (with "— >■ " indicating asymptotic values) 



c d 



e 



2(f> 



^ e-^^/c" (138) 



aciw 



w 



^ + e'^^) + hew 




(139) 



(140) 



f =p{c^w)+rc^{w'^ -e^'t') (141) 

^ 0. 

Here again, the inversion of e^^ indicates geometrically new AVTD solutions. 
As well, while the half-polarized solutions discussed above have w asymptot- 
ically vanishing and r approaching a general function on (see equations 
([57|)-(|58|)), these new solutions have r asymptotically vanishing and w ap- 
proaching a general function on T^. Examining the Hamiltonian (|])-(0), we 
see that this swapping of the asymptotic behavior of r and w is needed to 
avoid having H blow up as a consequence of the inversion of e^'^ . 



7 Conclusion 

How prevalent is AVTD behavior in solutions of Einstein's equations? Since 
this behavior is so specialized, it is perhaps rather surprising that, as we show 
here, it is found in substantial classes of solutions with only one Killing field. 
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As noted above, numerical work []^ indicates that in generic cosmological 
spacetimes with one Kilhng field, one finds Mixmaster rather than AVTD 
behavior near the singularity. However, it is likely that we can extend the 
classes of solutions known to have AVTD behavior in at least two ways. 

First, we should be able to remove the analytic condition which our the- 
orems here require. Rendall has shown how to do this for the class of Gowdy 
spacetime: In his work with Kichenassamy , it is shown that certain classes 
of analytic Gowdy spacetimes have AVTD behavior; then in his later work 
[llO| , he shows that the same holds for C°° solutions. Work hs begun which 
applies techniques similar to those used in |T0[ to show that AVTD behavior 
is found in classes of C°° polarized U{1) Symmetric spacetimes. 

Second, we hope to be able to show that all of the polarized U{1) Sym- 
metric vacuum solutions have AVTD behavior, and not just a subset of these 
spacetimes. Numerical studies suggest that this is true; we are looking for 
ways to prove this contention. 

It is not clear whether there are vacuum solutions with no Killing fields 
which exhibit AVTD behavior. Andersson and Rendall [|11|] have shown that 
there are classes of solutions of the Einstein equations coupled to a scalar 
field with AVTD behavior; however, while heuristic analyses have predicted 
this result, the same sorts of studies suggest that generic vacuum solutions 
with no Killing fields should show Mixmaster rather than AVTD behav- 
ior. We note that, to date, apart from the work of Berger and Moncrief 
PB| which uses Geroch transformations to generate a small class of U{\) 
Symmetric spacetimes with Mixmaster behavior, there are very few rigorous 
results concerning the existence of such behavior in spatially inhomogeneous 
cosmological spacetimes. 
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